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Abstract 


The Spike (S) protein of SARS-coronavirus (SARS-CoV) mediates viral entry into host cells. It contains two heptad repeat regions, 
denoted HRN and HRC. We have identified the location of the two interacting HR regions that form the six-helix bundle (B. Tripet, et al, 
J. Biol. Chem., 279: 20836-20849, 2004). In this study, HRC peptide (1150-1185) was chosen as the region to make structure-based substi- 
tutions to design a series of HRC analogs with increased hydrophobicity, helical propensity and electrostatic interactions, or with a cova- 
lent constraint (lactam bridge) to stabilize the a-helical conformation. Effects of the substitutions on o-helical structure of HRC peptides 
and their abilities to interact with HRN or HRC have been examined by biophysical techniques. Our results show that the binding of the 
HRC analogs to HRN does not correlate with the coiled-coil stability of the HRC analogs, but their interactions with HRC does corre- 
late with their stability, except for HRC7. This study also suggested three types of potential peptide inhibitors against viral entry can be 
designed, those that simultaneously inhibit interaction with HRC and HRN and those that are either HRC-specific or HRN-specific. For 
example, our study shows the important role of o-helical structure in the formation of the six-helix bundle where the lactam bridge con- 
strained analog (HRCS) provided the best interaction with HRN. The importance of o-helical structure in the interaction with native 
HRC was demonstrated with analog HRC4 which binds best to HRC. 
© 2006 Elsevier Inc. All rights reserved. 
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1. Introduction still remains a threat due to its highly transmittable nature 


to human populations and the mysterious origin of SARS- 


Severe acute respiratory syndrome (SARS) is an acute 
respiratory illness caused by infection with a novel corona- 
virus (SARS-CoV). The severity and mortality of this 
illness was witnessed during the global pandemic of SARS- 
CoV in 2003. The virus has not re-emerged since July 2003, 
except for several cases of laboratory-acquired infections 
and one natural outbreak resulting in four infected people 
(Normile, 2004a; Normile, 2004b). However, since SARS- 
CoV-related viruses have been detected in some animals, it 
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CoV (Guan etal. 2003; Hartley and Smith, 2003). Cur- 
rently, since there is no effective agent for the anti-viral 
therapy of SARS-CoV infection, it is imperative to learn as 
much as possible about this virus to accelerate the develop- 
ment of therapeutics and vaccines for its re-emergence. 
Infection by SARS coronavirus requires fusion of the 
viral and cellular membranes, which is mediated by 
the viral envelope Spike (S) glycoprotein and receptors on 
the target cell (Holmes, 2005). The S protein is a type I viral 
fusion protein which contains two highly conserved heptad 
repeat (HR) domains that have been shown to form coiled- 
coil structures (Bosch et al., 2003; Hakansson-McReynolds 
et al., 2006; Spaan et al., 1988; Tripet et al., 2004). The HR 
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region located closest to the N-terminus in the sequence is 
denoted HRN and the HR region located closest to the C- 
terminus is denoted HRC. Although the exact mechanism 
by which the SARS-CoV enters the host cell has not been 
elucidated, it is most likely similar to other coronaviruses. 
Upon binding to the receptor at the cell membrane, the 
fusion protein will be induced into the fusogenic intermedi- 
ate state with a dramatic conformation change. Collapse of 
this fusogenic intermediate leads to a six-helix bundle (tri- 
mer of dimers) formation between HRN and HRC, which 
ultimately promotes membrane fusion (Eckert and Kim, 
2001b). Our laboratory and other researchers recently 
reported (Bosch et al., 2004; Ingallinella et al., 2004; Tripet 
et al., 2004; Xu et al., 2004b) that HR regions of SARS- 
CoV S protein can associate to form a very stable o-helical 
six-stranded structure and the orientation of the HR 
regions is anti-parallel. Residues 902-950 in the HRN 
region and 1151-1185 in the HRC region were identified to 
be crucial for their interaction. The structure of the ectodo- 
main of SARS Spike protein in its fusogenic/post-fusogenic 
state and HR interaction regions have recently been con- 
firmed by crystallography (Duquerroy et al., 2005; Supekar 
et al., 2004; Xu et al., 2004a). Most recently, Li et al. (2005) 
reported the crystal structure of the SARS-CoV S receptor 
binding domain (residues 306-575) complexed with recep- 
tor ACE2 (angiotensin-converting enzyme 2), which further 
confirms the prediction that SARS-CoV fusion process is 
similar to other viruses with type I fusion proteins. Hakans- 
son-McReynolds et al. (2006) also reported the solution 
structure of the SARS-Coronavirus HRC domain (1143— 
1193) in the prefusion state, which forms a coiled-coil sym- 
metric trimer. 

A critical step in viral entry is the conformational change 
of the fusion domain of the S protein. Peptides derived 
from the HRC region of the HIV-1 gp41 and other class I 
fusion proteins have been reported to show significant viral 
fusion inhibitory activity (Chan and Kim, 1998; Eckert and 
Kim, 2001a; Root et al., 2001). It has been proposed that 
HRC analogs bind strongly to the transiently exposed 
HRN coiled-coil trimer, and block the formation of the six- 
helical bundle necessary for the fusion. For HIV and 
murine CoV mouse hepatitis virus (MHV), HRC peptides 
have been shown to inhibit viral entry at nanomolar con- 
centrations (Bosch etal. 2004; Chan and Kim, 1998). 
Recently, a HIV peptide inhibitor has been made by the 
fusion of gp41 HRN region and a trimeric forming coiled- 
coil peptide, where the three chains are covalently linked by 
disulfide bridges. This construct was a potent inhibitor of 
viral entry at pM concentration (Bianchi et al., 2005). One 
of the HRC peptides of HIV-1 gp41, Enfuvirtide (Fuzeon), 
was approved by the FDA for treatment of AIDS (LaBo- 
nte et al., 2003). In several recent studies (Bosch et al., 2004; 
Liu et al., 2004; Yuan et al., 2004; Zhu et al., 2004), SARS- 
CoV S-mediated fusion was shown to be inhibited by 
HRC-derived peptides; however, these peptides inhibited 
viral entry at concentrations in the micromolar range 
(Bosch et al., 2004; Liu et al., 2004). One possible reason for 


the higher inhibitory concentration is that the interaction 
between the inhibitory peptides and HRN or HRC is too 
weak to prevent the native intra-molecular association of 
HRN and HRC. As reported in studies of HIV peptide 
inhibitors, considerable genetic and biophysical evidence 
supports the concept that the ability of the Class I envelope 
glycoprotein to mediate membrane fusion is determined, in 
part, by the high stability of the six-helix bundle (Eckert 
and Kim, 2001b; Skehel and Wiley, 2000). Thus, peptides 
inhibiting the fusion process must be designed with the 
highest binding affinity/stability to be useful for preventing 
viral entry. 

In the present study, HRC peptide (1150-1185) was cho- 
sen as the region to make structure-based substitutions to 
design a series of HRC analogs. These HRC analogs were 
designed with increased hydrophobicity, helical propensity 
and electrostatic interactions, or with a covalent constraint 
(lactam bridge) to stabilize the o-helical structure. Effects of 
the substitutions on the conformation of HRC peptides 
and their abilities to interact with HRN or HRC have been 
examined by biophysical techniques. 


2. Materials and methods 
2.1. Preparation of peptides 


2.1.1, Peptide synthesis 

The HRN and HRC peptides of SARS-CoVS glycopro- 
tein were prepared by solid-phase synthesis methodology- 
using 4-benzylhydrylamine hydrochloride resin with 
conventional fmoc (fluorenylmethoxycarbonyl) chemistry 
as described by Tripet et al. (2000). 

The synthesis of the covalent HRC trimer and HRN tri- 
mer for surface plasmon resonance analysis was carried out 
as described by Tripet et al. (2006). 

The gene sequence Gene bank Accession No. for the 
SARS-CoV S protein is AY 278741. 


2.1.2. Lactam formation on peptide resin 

The protecting groups for the side-chains of Lys1170 
and Glu1166 in the HRCS sequence which forms the 7, i+ 4 
lactam bridge were the allyloxycarbonyl group (Aloc) for 
Lys and the allylester for Glu, Glu (OAI). These two 
protecting groups were removed by Tetrakis(triphenyl- 
phosphine)palladium(0) [Pd(PPh3),] as described by Kates 
et al. (1993). The reaction was carried our under argon. 
Formation of the lactam ring was achieved with benzotriaz- 
ole-1-yl-oxy-tris(pyrrolidino)phosphonium_ hexafluorophos- 
phate (pyBOP) in the presence of N,N-diisopropylethylamine 
(DIEA). 


2.1.3. Peptide purification and characterization 

Peptides were N-terminally acetylated, cleaved from the- 
resin, and purified by reversed-phase high-performance liq- 
uid chromatography (RP-HPLC) on a Zorbax 300 SB-C8 
preparative column (250mm x 9.4mm LD., 6.5 um particle 
size, 300-A pore size; Agilent Technologies, Little Falls, 
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DE) with a linear AB gradient (0.6% B/min) at a flow rate 
of 2 ml/min, where eluent A is aqueous 0.05% trifluoroace- 
tic acid (TFA) and eluent B is 0.05% TFA in acetonitrile. 
Homogeneity of the purified peptides was characterized by 
analytical RP-HPLC, amino acid analysis, and electrospray 
mass spectrometry. 


2.2. Circular dichroism spectroscopy 


Circular dichroism spectra were recorded on a Jasco 
J-810 spectropolarimeter (Jasco Inc., Easton, MD). The CD 
wavelength scans were measured from 190 to 255nm in 
benign buffer (0.1 M KCl, 0.05 M PO,, pH 7.0). For samples 
containing trifluoroethanol (TFE), the above buffer was 
diluted1:1 (v/v) with TFE. Temperature denaturation mid- 
points (7,,) for the peptides were determined by following 
the change in molar ellipticity at 222 nm from 4 to 95°C in 
a 5-mm path length cell and a temperature increase rate of 
1°C/min. Ellipticity readings were normalized to the frac- 
tion of peptide folded (f) or unfolded (f,) using the stan- 
dard equations f-= ((0] — [6],,)/({@], —[6],) and f,=( —/fp; 
where [6],, and [0], represent the molar ellipticity values for 
the native fully folded and fully unfolded species, respec- 
tively. [6] is the observed molar ellipticity at 222 nm at any 
temperature. 


2.3. Size-exclusion and reversed-phase HPLC 


HRN and HRC peptides singly or as an equimolar mix- 
ture (200M of each peptide) were dissolved in 100mM 
KCl, 50mM PO, buffer, pH 7, and equilibrated at room 
temperature for 30min. A 10-ul aliquot of sample was ana- 
lyzed with a high-performance size-exclusion column, 
Superdex 75™ (30cm x 1cm I.D., Amersham Biosciences, 
Uppsala, Sweden) equilibrated in a buffer consisting of 
50mM sodium phosphate/150mM NaCl, pH 7.5, at a flow 
rate of 0.5 ml/min and ambient temperature. The molecular 
weights of the peak fractions were estimated by comparison 
with a protein standard running on the same column. 

For peptides that formed stable hetero-stranded com- 
plexes, the complex peak was collected and analyzed by 
RP-HPLC on an analytical C8 column (Zorbax 300SB-C8, 
150mm x 4.6mm LD, 6.5 um particle size, 300-A pore size). 
The peptides were eluted from the column by employing a 
linear AB gradient of 2% B/min, where eluent A is 0.05% 
aqueous TFA and eluent B is 0.05% TFA in acetonitrile at 
a flow rate of 1.0 ml/min at room temperature. 


2.4. Native-polyacrylamide gel-electrophoresis (PAGE) 


Equimolar mixtures (200 uM of each peptide) of HRN 
(in phosphate buffer, pH 7.3) and each HRC analog were 
incubated at 25°C for 10min (sample volume was 51). 
After the addition of an equal volume of 2 x native sample 
buffer (0.125M Tris-HCl, pH 6.8/10% glycerol/0.004g of 
bromophenol blue) the mixture was analyzed by PAGE on 
a 15% Tricine gel with a Tricine/glycine running buffer 


(pH 8.3). Gel electrophoresis was carried out with 120V 
constant voltage at 4°C for 2h. The gel was then stained 
with Coomassie blue. 


2.5. Surface plasmon resonance (SPR) analysis 


To measure the binding affinity or kinetic dissociation 
constant of the SARS-CoV related peptides, we used sur- 
face plasmon resonance with an optical biosensor (model 
BIACORE 3000, BIACORE Inc, Piscataway, NJ, USA). 

Covalent coupling of HRN trimer or HRC trimer to the 
BIACORE CMS chip (a dextran-coated sensor chip) sur- 
face was performed following the standard amine coupling 
procedure according to the manufacturer’s specification. 
Briefly, the dextran surface of the CM5 chip was first 
activated with 20 yl N-hydroxysuccinimide (NHS)/N-ethyl- 
N'-(dimethyl-aminopropyl)-carbodiimide (EDC). HRN tri- 
mer or HRC trimer in 10mM sodium acetate buffer 
(pH4.3) was injected and allowed to react to give a surface 
density of 420 resonance units (RU) for the HRC trimer 
chip and 2000 RU for the HRN chip. Remaining activated 
groups were then blocked by injection of 35 ul ethanol- 
amine HCl. A reference surface with the same treatment 
except no injection of ligand was generated at the same 
time and was used as background to correct instrument and 
buffer artifacts. 

Binding experiments were performed at 25°C in PBS 

buffer, pH 7.4, with 0.005% Tween 20. The peptides to be 
tested were introduced at four serial concentrations (5, 10, 
20, and 50 1M) in running buffer in 80 ul injection volumes 
and at a flow rate of 20 l/min, followed by 6 min of dissoci- 
ation. Each experiment was carried out at least three times. 
Remaining analytes were removed in the surface regenera- 
tion step using 6M guanidine hydrochloride in 20mM gly- 
cine, pH 2.0 for the HRN trimer chip, and 100mM acetate 
buffer, pH 4.0 for the HRC trimer chip. 
In Tables 2 and 4, the K,,, and Ko values could not deter- 
mined since the simple 1:1 binding model could not be used 
to fit the observed data. The results are complicated by the 
fact that the HRC analog flowing over the chip surface is in a 
monomer-oligomer equilibrium and the binding event on the 
immobilized HRC trimer or HRN trimer could involve the 
binding of one or more HRC polypeptide chains. At present 
no model is available to calculate these parameters for such a 
complex system. Thus, a simple “+” or “—” was used to 
denote binding or no binding, respectively. 


3. Results 
3.1. Peptide design 


Residues 902-950 in the HRN region and 1150-1185 in 
the HRC region were identified to be crucial for their inter- 
action to form the six-helix bundle state (Bosch et al., 2004; 
Ingallinella et al., 2004; Tripet et al., 2004; Xu et al., 2004b). 
Furthermore, it has been reported that SARS-CoV S-medi- 
ated fusion can be inhibited in a concentration-dependent 
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manner by HRC-derived peptides in this HRC region (Liu 
et al., 2004; Zhu et al., 2004). Therefore, in our study, the 
36-residue HRC peptide (1150-1185) was chosen as the 
region to make structure-based substitutions to design a 
series of HRC analogs, to increase coiled-coil stability and 
their ability to interact with HRN (902-950) or HRC 
(1150-1185) (Fig. 1). 

In the design of peptide inhibitors, one must consider the 
following: first, in the current model of membrane fusion, 
the anti-viral activity of peptide fusion inhibitors would 
depend on the binding affinity between the peptide and the 
HRN region or HRC region of SARS-CoV S82. Second, 
substitutions in the peptide sequence to change its confor- 
mation or stability could substantially affect the interaction 
with the target (HRN or HRC) (Chan et al., 1998; Eckert 
and Kim, 2001a; Lu et al., 1995). For example, the HRC 
region was shown to have o-helical structure on its own 
(Hakansson-McReynolds et al., 2006; Tripet et al., 2004), 
suggesting that enhancement of o-helical structure and sta- 
bility could increase affinity of potential peptide inhibitors 
for their helical targets. Factors that contribute to the heli- 
cal conformation of peptides and stability of coiled-coils 
have been systematically studied (Chakrabartty et al., 1994; 
Houston et al., 1995; Houston et al., 1996; O’Neil and DeG- 
rado, 1990; Tripet et al., 2000; Wagschal et al., 1999b; Zhou 
et al., 1994). There are two general ways to stabilize the 
coiled-coil. Increasing the hydrophobicity of the residues in 
the hydrophobic core stabilizes the coiled-coil (Tripet et al., 
2000; Wagschal et al., 1999b) or increasing interactions that 
stabilize the individual o-helix (increasing o-helical propen- 
sity of residues, adding intra-chain electrostatics or intra- 
chain covalent constraints) which in turn stabilizes the 
coiled-coil by shifting the monomer-trimer equilibrium in 
favor of the folded trimeric coiled-coil. Thus, we designed 
the substitutions shown in Fig. 1 which involved: 


(1) Increasing hydrophobicity in the hydrophobic core: 
(HRC1). 
(2) Increasing helical propensity (HRC2 and HRC4). 
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(3) Increasing hydrophobicity in the hydrophobic core 
and increasing helical propensity (HRC3). 

(4) Introducing an i to i+ 4 lactam bridge (HRC%S). 

(5) Introducing intrachain i to i+3 or i to i+4 salt 
bridges (HRC6 and HRC7). 


3.2. Circular dichroism analysis of HRC analogs 


To examine the change in o-helical structure of the HRC 
analogs, we analyzed each peptide by circular dichroism 
spectroscopy (CD). As seen in Fig. 3, under physiological 
conditions, all of these analogs are well folded, soluble, and 
have typical o-helical spectra. 

The content of a-helicity of the native HRC peptide 
was concentration-dependent. It appeared to be largely 
unfolded at low concentrations (26% helix at 7uM), 
while at 70 uM it was 65% o-helical (Fig. 2 and Table 1). 
It was previously shown by sedimentation equilibrium 
experiments that HRC forms a trimer in the analytical 
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Fig. 2. CD spectra of native HRC (1150-1185) at different concentrations. 
Spectra were recorded in a 0.1 M KCI, 0.05 M PO, buffer, pH 7. 


902 950 
d ada ad adaadaal| 
HRN Ac-QKQIANOFNKAISQIQESLTTTSTALGKLODVVNONAQALNTLVKQLSS-amide 
. = . = . 
Region1 +! Region 2 : Region 3 
1150 - 1185 
ad ae! a aa ad 
HRC Ac-DISGINASVVN » IQKEIDRLNEVAKNL = NESLIDLQEL-amide 
HRC1 Ac-DISGINASVVN § IQKEIDRLNEVIKNL § NESLIDLQEL-amide 
HRC2 Ac-DISGINASVVN § IOQKEIARLNEVAKAL © NESLIDLQEL-amide 
HRC3 Ac-DISGINASVVN ® IQKEIARLNEVIKAL * NESLIDLQEL-amide 
HRC4 Ac-DIAAINASVAN # IQKEIARLNEVAKAL # NESLAALQAL-amide 
1 1 
HRC5 Ac-DISGINASVVN =» IQKEIERLNKVAKNL = NESLIDLQEL-amide 
HRC6 Ac-DISGINASVVN « IQKEIERLNKVAKNL » NESLIDLQEL-amide 
HRC7 Ac-DISGINASVVE § IOKKIEELNKKAEEL 3 NKKLIDLQEL-amide 
a 


Fig. 1. Amino acid sequences of peptides used in the study. The a and d positions of the predicted coiled-coil heptad repeats (abcdefg)n are shown above 
the sequence. HRCI—HRC7 are HRC-derived peptides, where the residues in bold and underlined are substitutions in the native HRC (1150-1185) 
sequence. The HRC sequence was divided into 3 regions based on the three-dimensional structure of HRC in the six-helix bundle state, Region | and 3 are 


extended in the structure and Region 2 is o-helical. 
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Table 1 
Ellipticity and stability of the synthetic peptides 


Peptide name [O]x09" % o-helix® [0}222/208 Tm (°C) Oligomerization state‘ 
Benign 50% TFE Benign 50% TFE Benign 50% TFE 

HRC —21370 —29720 65 91 0.86 0.85 37 Trimer 

HRC 1 —32060 —31510 98 96 0.97 0.86 57 Tetramer 

HRC 2 —27060 —29 740 83 91 0.93 0.85 48 Trimer 

HRC 3 —30900 —31 450 95 96 0.98 0.86 66 Tetramer 

HRC 4 —28 560 —30 990 87 95 1.09 0.90 74 Tetramer 

HRC 5 —27770 —32 320 85 99 0.87 0.86 41 Trimer 

HRC 6 —16510 —30 190 51 92 0.79 0.84 32 Trimer 

HRC7 —32440 —32 660 99 100 1.01 0.88 61 Trimer 


4 [0]. is the mean residue molar ellipticity (degrees cm? dmol™!) measured at 222 nm in a 0.1 M KCl, 0.05 M PO, buffer, pH 7, in the absence (Benign) or 
presence of 50% trifluoroethanol (50% TFE) (v/v). Concentration of peptides was 70 1M. 

> The ellipticity value in 50% TFE for peptide HRC7 was taken as 100% and % a-helical content of the other peptides were relative to HRC7. 

© T,, is the temperature at which there is a 50% decrease in o-helical content compared to the fully folded coiled-coil as determined by circular dichroism 


spectroscopy at 5 °C. 


4 Oligomerization state of peptides was determined by SEC-HPLC. Concentration of peptides was 100 uM. 


ultracentrifuge (Tripet et al., 2004). The CD results agree 
with a monomer-trimer equilibrium at low concentra- 
tions and a fully-folded trimer state occurring at 70 uM. 
Thus, the HRC peptide does not form a very stable 
oligomer. 

As illustrated in Fig.3 and Table 1, the HRC analogs 
had a higher content of helical structure compared to HRC, 
except for HRC6 which was designed as the control peptide 
for HRC5 and HRC7, showed lower helical structure than 
native HRC. HRC1, HRC3, and HRC7 are essentially fully 
folded in benign buffer (>95% a-helical content, Table 1) 
and thus no significant change in helical structure is 
observed by addition of the helix inducing solvent, TFE. 
However, HRC and HRC6 were not fully folded in benign 
buffer and substantial o-helical structure could be induced 
by addition of TFE. 

To test the stability of the HRC analogs, each analog 
was denatured by increasing temperature, as illustrated in 
Fig. 4, and their 7,, values are shown in Table 1. These 
results showed that all analogs were more stable than 
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Fig. 3. CD analysis of helical structure of HRC analogs at 70 uM. Spectra 
were recorded in a 0.1 M KCl, 0.05 M PO, buffer, pH 7. 
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Fig. 4. Stability of HRC analogs. Temperature denaturation profiles of the 
HRC analogs monitored by CD at 222 nm in a 0.1M KCl, 0.05M PO, 
buffer, pH 7. Concentration was 701M. The fraction folded (f;) of each 
peptide was calculated as f;= ((0]—[0],)/((@], —[@],), where [0] is the 
observed mean residue molar ellipticity at 222 nm at any particular tem- 
perature, and [6], and [6], are the mean residue molar ellipticities at 
222 nm of the native folded state at 4°C and unfolded states, respectively. 


native HRC, as expected, except for HRC6, which was less 
stable than native HRC (7,,, values range from 32 to 74°C). 

Overall, the change in structure of HRC with substitu- 
tions was consistent with our design aim. That is the stabil- 
ity of the HRC1 coiled-coil increased due to the increase in 
hydrophobicity in the hydrophobic core by the Ala to Ile 
substitution at position 1172a (Fig. 1), as suggested by the 
hydrophobicity coefficients of all 20 amino acid residues in 
positions a and d of coiled-coils determined by Hodges and 
coworkers (Tripet et al., 2000; Wagschal et al., 1999a). Sub- 
stitutions made to residues occurring in the b, c, and f hep- 
tad positions (the outer surface of the coiled-coil) with Ala, 
which has the highest a-helical propensity of all side-chains 
(Zhou et al., 1994), increased the a-helical structure and 
stability of HRC2 and HRC4 coiled-coils (Figs. 3 and 4, 
Table 1). The larger number of Ala substitutions in HRC4 
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substantially increased the stability relative to HRC2 (8 vs 
2 Ala substitutions resulted in T,, values of 74 and 48°C, 
respectively; Table 1). Peptide HRC3, which contained 
both the Ie residue in the hydrophobic core and the two 
Ala residues formed a more stable coiled-coil than either 
HRCI and HRC2 analogs. A comparison of the substitu- 
tions showed that Ile resulted in a greater increase in coiled- 
coil stability than the two substituted Ala residues (AT,, 
HRC1-HRC was 21°C compared to AT,, HRC2-HRC of 
11°C) indicating, as expected, that a single substitution in 
the hydrophobic core made a greater contribution to 
coiled-coil stability compared to two Ala substitutions in 
positions on the outer surface of the coiled-coil. However, 
the a-helical propensity effect for increasing stability can 
dominate over the increase in stability observed for the sin- 
gle Ala to Ie substitution in the hydrophobic core if the 
number of substitutions increasing a-helical propensity is 
large enough. That is, the 8 Ala substitutions (HRC4) 
increased stability by 37°C (AT,, HRC4—HRC) compared 
to the single Ile substitution in the hydrophobic core (AT, 
HRC1-HRC) of 21°C. 

The formation of the lactam bridge on the outer face of 
each helix significantly stabilized the coiled-coil of HRC5 
compared to the non-covalent salt bridge in HRC6 (7, of 
32°C for HRC6 compared to a T,, of 41°C for HRC5; 
Fig. 4 and Table 1). Intrachain electrostatics also modu- 
lated stability of helical structure. Substitutions in HRC6 
introduced an intrachain i to i+ 4 ionic attraction between 
E1156 and K1160 that was expected to increase o-helical 
structure and stability. However, K1160 also introduces 
two potential i to i+3 charge repulsions with R1157 and 
K1163 which could dominate over the introduced attrac- 
tion. Thus, HRC6 which was destabilized relative to HRC 
(AT, of 5°C) was a control peptide for the covalent lactam 
bridge peptide HRCS. When a series of intrachain i to i+3 
and i to i+4 salt bridges were introduced into HRC, the 
resulting analog HRC7 (Fig. 1) had a stability of 61°C 


compared to HRC6 with the single salt bridge addition 
(T,, 32°C) or HRC (T,, 37°C) (Table 1). 

The oligomerization state of these analogs was deter- 
mined by size-exclusion chromatography (SEC) and shown 
to be either trimers or tetramers (Table 1). The oligomeriza- 
tion specificity of coiled-coils can be switched by even a sin- 
gle substitution, e.g., HRC (trimer) is switched to HRC1 
(tetramer) by A1172I substitution. 


3.3. Interaction of HRC analogs with HRN 


3.3.1. Circular dichroism analysis 

We used CD to analyze the ability of these HRC analogs 
to bind HRN. The results (Table 2) showed that interaction 
with HRN induces a significant change in o-helical confor- 
mation (e.g., HRN + HRCS, shown in Fig. 5), as compared 
with a theoretical spectrum obtained by summing experi- 
mental spectra of equivalent amounts of HRN and HRC 
analogs alone. These results suggested that HRC and ana- 
logs HRC1, 2, 3, 5, and 6 interact with HRN, whereas 
HRC4 and 7 have weak or no interaction (Table 2). 

To test the stability of complexes of HRN with HRC 
analogs, 1:1 molar equivalent complexes of HRN and HRC 
or HRC analogs were determined by increasing tempera- 
ture, as illustrated in Fig. 6. The results show that HRC 
(native sequence) can form a very stable complex with 
HRN, which agrees with our choice of the HRC sequence 
(1150-1185) as the key region of interaction with HRN and 
is suitable for substitutions to design analogs with 
improved properties. 

A shown in Table 2, the substitutions resulted in a set of 
HRC analogs that form a set of HRC/HRN complexes 
with 7, values ranging from 59 to 88 °C compared to HRN 
alone with a T,, of 43 °C, except for HRC7 where the T,, of 
the mixture with HRN was similar to HRN or HRC7 
alone, indicating that HRC7 does not interact with HRN 
(Fig. 6). Native HRC interacted and formed a more stable 


Table 2 

Summary of interaction of HRC analogs with HRN 

Complex [O]229 Change in [0] T° (°C) AT,,° Complex formation tested by SPR binding 

e 

Theoretical* Observed HRN HRN-HRC_ Native-PAGE® SEC-HPLC? *st 

HRN —20 890 43 

HRN-—HRC —15920 —19060 3140 85 +42 0 + + + 

HRN-HRC1 —23790 —25 340 1540 74 +31 0-11 + + + 

HRN-HRC2 —19060 —23070 4010 86 +43 +1 + + + 

HRN-—HRC3 = —24 250 —26 130 1880 75 +32 —10 + + + 

HRN-HRC4 —25610 —26 180 570 59 +16 —26 ~ + + 

HRN-HRC5 = —22420 —27700 5280 88 +45 +3 + + + 

HRN-HRC6 —15980 —21180 5200 83 +40 —2 + + + 

HRN-HRC7  —23 340 —23 850 510 44 +1 —41 — — 


* The theoretical molar ellipticity for two non-interacting peptides was calculated by summing the two individual spectra. 

2 T, 1s the transition midpoint temperature at which there is a 50% decrease in molar ellipticity [0], compared with the fully folded peptide as deter- 
mined by CD at 5 °C. The CD test was carried out after incubating the peptide mixture for 20 min at 5 °C. Concentration of each peptide was 7 1M. 

° AT,, is the change in T,, from the T,, of HRN or the T,, of HRN—HRC complex. 

4 HRC peptides were incubated with HRN for 20 min, and then analyzed by Native-PAGE and SEC-HPLC. Concentration of each peptide in the mix- 


ture was 100 uM. 


° SPR binding analysis was carried out on the CMS sensor chip with the immobilization of the HRN trimer peptide (see Section 2). 
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Fig. 5. CD analysis of the interaction of HRN and HRCS. CD spectrum 
of a 1:1 molar complex between HRN and HRC5 at 25 °C ina 0.1 M KCl, 
0.05 M PO, buffer, pH 7. Peptide concentrations were 7 uM. The theoreti- 
cal spectrum for two non-interacting peptides is shown for comparison, 
which was generated by adding the individual peptide spectra at the same 
concentrations. 


complex with HRN than peptide analogs HRCI, HRC3, 
and HRC4. However, analogs HRC2, HRCS5, and HRC6 
formed complexes as stable as or slightly more stable than 
native HRC/HRN. HRN/HRC5 complex was not com- 
pletely unfolded, even at 95°C (Fig. 6). 


3.3.2. Size-exclusion chromatography (SEC) analysis 

To show complex formation between HRN and HRC 
analogs, we collected the earlier eluted peak from SEC and 
resolved it by RP-HPLC (Fig. 7). As shown in Table 2, 
except for HRC7, the other HRC analogs form stable com- 
plexes with HRN. The elution volumes of all the complexes, 
relative to Mr standards, gave an apparent M, of 27490 Da, 
in good agreement with the M, (27 900-28 340 Da) expected 
for a hexameric complex, formed by three molecules of 
HRN and three molecules of HRC or HRC analogs. 
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Fig. 7. SEC analysis of the HRN and HRC complex. HRN (2 nmol) and 
HRC (2 nmol) were pre-incubated together for 30 min in running buffer 
and then 10 tl of incubated sample was applied to SEC; conditions shown 
in Section 2.3. The absorbance peak at 19.7 min corresponding to the 
HRN/HRC complex was collected and subsequently analyzed by RP- 
HPLC (inset); Conditions shown in Section 2.3. Each absorbance peak is 
labeled accordingly. 


As shown in Fig. 7 (inset), the integrated peak areas con- 
verted to concentration using standard curves showed that 
the two peaks of HRN and HRC were in a 1:1 mole ratio as 
determined by RP-HPLC. 


3.3.3. Native-polyacrylamide gel-electrophoresis (PAGE) 
analysis 

To verify further complex formation between HRN and 
HRC analogs, we analyzed their interaction by Native- 
PAGE. The results are shown in Fig. 8 and Table 2. HRC 
or HRC analogs alone showed a band in the lower part of 
the gel. With the exception of HRC4 and HRC7, the mix- 
ture of HRC analogs and HRN showed two bands: the 
lower band had the same position as isolated HRC analog, 
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Fig. 6. Temperature denaturation profiles of HRN alone and a 1:1 molar mixture of HRN with HRC analog monitored by CD at 222 nm ina 0.1 M KCl, 
0.05 M PO, buffer, pH 7. Concentrations of peptides were 7 1M. Fraction folded was calculated as described in Fig. 4 legend. 
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Fig. 8. Complex formation between HRC analogs and HRN, as deter- 
mined by Native-PAGE. HRN and HRC analogs on their own or as pre- 
incubated equimolar (200 4M of each peptide) mixture were subjected to 
15% PAGE. Samples were incubated for 30 min in 0.1 M KCl, 0.05 M PO, 
buffer, pH 7, diluted 1:1 (v/v) with native sample buffer and loaded onto 
the gel. 


and the upper band is higher order oligomeric complex 
formed by HRN and HRC analog. HRN alone showed no 
band because it carries a net positive charge under the 
native gel electrophoresis conditions, and consequently 
does not enter the gel. 


3.3.4. Surface plasmon resonance analysis of HRC analogs 
binding HRN 

We next carried out SPR binding analysis of these HRC 
analogs on BIACORE chips with the immobilization of 
HRN trimer peptide. As shown in Fig. 9 and Table 2, with 
the exception of HRC7, HRC, and HRC analogs could 
bind the HRN trimer chip, which is consistent with the 
results from SEC analysis. The regeneration step for the 
chip surface required the use of a strong denaturant (6M 
guanidine hydrochloride in 20mM glycine, pH 2.0 for 
1 min). This treatment also showed HRC analogs bind the 
HRN trimer very strongly. 
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Based on the results in Table 2 obtained from CD, SEC, 
Native-PAGE, and SPR analysis, compared to native 
HRC, the following results were obtained: HRC2 and 
HRCS5 have high affinity for HRN; HRC1 and HRC3 
showed decreased interaction with HRN and HRC4 could 
only form a weak complex with HRN. All the data showed 
no interaction between HRC7 and HRN. The results illus- 
trate that when substitutions increase the stability of HRC 
helical structure (Table 1), such as that seen for HRCI, 
HRC3, and HRC4 (Table 2), the stability of the complex 
with HRN decreases. Therefore, the binding affinity of the 
analogs for HRN does not show a clear correlation to the 
coiled-coil stability of the HRC analogs. Even though HRC 
and HRC6 (Table 1) are only partially folded o-helical 
structures at the concentration of 71M, they can form very 
stable complexes with HRN (7, of 85 and 83°C, respec- 
tively, Table 2). 

In this study, the HRC sequence was shortened to 28 res- 
idues (1158-1185) and 25 residues (1161-1185). The results 
in Table 3 showed that the 28-residue peptide had a much 
weaker interaction with HRN (a complex with a T,, of 
58°C) than the 36-residue HRC (1150-1185) (a complex 
with a T., of 85°C). When shortened to 25 residues (1161— 
1185), this analog showed no interaction with HRN. Even 
after formation of one lactam bridge was made in the core 
of the 25-residue HRC peptide to stabilize a-helical struc- 
ture, it still did not show any interaction with HRN. This 
result is in good agreement with the previous reports 
(Bosch et al., 2004; Ingallinella et al., 2004; Zhu et al., 2004). 


3.4. Interaction of HRC analogs with HRC 


3.4.1. CD analysis 

As shown in Fig.4 and Table 1, the stability of most 
HRC analogs increased in comparison to the native HRC. 
One can propose that even though these HRC analogs are 
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Fig. 9. Surface plasmon resonance analysis of the interaction of HRC6 and HRC7 with a HRN trimer chip. HRN trimer was covalently coupled to the 
Biacore CMS chip surface. Binding experiments were performed at 25°C in PBS buffer, pH 7.4, with 0.005% Tween 20. The peptides to be tested were 
introduced at four serial concentrations (5, 10, 20, and 50 1M) in running buffer in 80 pl injection volumes and at a flow rate of 20 l/min, followed by 
6 min of dissociation. The schematic shows HRC6 or HRC7 as a monomer-trimer equilibrium in solution when flowing over the immobilized HRN trimer 


on the chip surface. 
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Table 3 
The structure of short HRC analogs and their interaction with HRN? 


Short HRC analogs Sequence a-Helical structure = Interaction with HRN 

CD Native-PAGE SEC _ 7, (°C) 
HRC 25 Ac-IQKETDRLNEVAKNLNESLIDLQEL-amide No a = a = 
HRC 28 Ac-VVNIQKEIDRLNEVAKNLNESLIDLQEL-amide No + + + 58°C 
HRCS 25 Yes - — - - 


Ac-IQKEIERLNKVAKNLNESLIDLQEL-amide 


@ The conditions for CD, Native-PAGE, and SEC were the same as shown in Table 2. 


more stable than HRC, they could still interact with native 
HRC to form mixed oligomers and thus prevent the forma- 
tion of the hexameric complex between native HRC and 
HRN. To verify this assumption, we first used CD to ana- 
lyze the interaction of these HRC analogs with native 
HRC. The results in Table 4 showed that the change in heli- 
cal structure is induced by the interaction of HRC analogs 


with HRC, which suggested that HRC analogs HRCI, 2, 3, 
4, and 5 interact with native HRC, whereas HRC6 and 7 
have weak or no interaction. 

To test if HRC analogs can stabilize HRC structure, 
temperature denaturation profiles were determined for a 
1:1 mixture of HRC and HRC analogs. The results are 
illustrated in Fig. 10, and 7,, values are shown in Table 4. 


Table 4 

Interaction of HRC analogs with HRC 

Complex name [O29 Change in [6])95 Tm? (°C) Interaction tested SPR binding test‘ 

ai ic 

Theoretical* Observed py Neve eee 

HRC 12320 5 

HRC-—HRC1 —19650 —25 240 5590 41 ND + 

HRC-—HRC2 —11090 —16280 5190 34 + + 

HRC-—HRC3 —19000 —26840 7840 46 + + 

HRC—HRC4 —20 600 —25 120 4520 55 + + 

HRC-—HRCS5 —18470 —21790 3320 38 + + 

HRC—HRC6 —9240 —10610 1370 5 — + 

HRC—HRC7 —18 700 —20 280 1580 42 - _ 


* The predicted molar ellipticity for two non-interacting peptides was calculated by summing the two individual spectra. 

b T, is the transition midpoint temperature at which there is a 50% decrease in molar ellipticity [0]. compared with the fully folded peptide as deter- 
mined by CD. Concentration for each peptide was 71M. The CD determination was carried out after incubating the peptide mixture for 20 min at 4 °C. 

° HRC analogs were incubated with native HRC for 20 min, and then analyzed by Native-PAGE. Concentration of each peptide for forming complexes 
was 100 uM. ND denotes that the interaction between HRC and HRCI could not be determined by this method since these two peptides show the same 


mobility in gel-electrophoresis. 


4 SPR binding analysis was carried out on a CMS chip with the immobilization of HRC trimer peptide (see Section 2). 
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Fig. 10. Temperature denaturation profiles of HRC alone and a 1:1 molar mixture of HRC with HRC analog monitored by CD at 222 nm in a 0.1 M KCl, 
0.05 M K,HPO, buffer, pH 7. Peptide concentrations were 7 .M. Fraction folded was calculated as described in Fig. 4 legend. 
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These results show that HRC3 and HRC4 had the stron- 
gest interaction with HRC. Thus, the substitutions resulted 
in a set of HRC analogs that interact with native HRC with 
a set of 7,,, values ranging from 34 to 55°C compared to the 
stability of HRC alone with a T,, of 5°C. The T,, for the 
HRC6 and HRC mixture is similar to native HRC, suggest- 
ing no interaction with HRC. 


3.4.2. Native-PAGE analysis 

In the gel-shift experiment, the complex formation of 
HRC with HRC analogs was observed on native gel-elec- 
trophoresis to compare with the HRC and HRC analogs 
alone as shown in Fig. 11 and Table 4. HRC or HRC ana- 
logs alone showed a single band. With the exception of 
HRCI, HRC6, and HRC7, the mixture of HRC analogs 
and HRC showed smeared bands between the positions of 
the two bands of native HRC and HRC analog. These 
smeared bands may be caused by the equilibrium between 
bound and unbound and various combinations of HRC 
and HRC analog (3:0, 1:2, 2:1, and 0:3, respectively) 
would have different mobilities due to different mass/ 
charge ratio. HRC] and HRC showed the same mobility 
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Fig. 11. Native-PAGE analysis of interaction of HRC analogs with HRC. 
HRC and HRC analogs on their own or as pre-incubated equimolar 
(200 uM of each peptide) mixture were subjected to 15% PAGE. Samples 
were incubated for 30 min in 0.1 M KCl, 0.05 M PO, buffer, pH 7, and the 
diluted 1:1 (v/v) with native sample buffer and loaded onto the gel. 


Response 


Time (sec) 


in the gel; thus, the interaction between them could not be 
determined by this method. However, the mixture of HRC 
and HRC6 or HRC7 showed two single bands corre- 
sponding to HRC and HRC6 or HRC7 alone, which sug- 
gests no interaction between HRC and HRC6 or HRC 
and HRC7 (Fig. 11). 


3.4.3. Surface plasmon resonance analysis of HRC analogs 
binding HRC 

We have carried out SPR binding analysis of these HRC 
analogs on BIACORE chip with the immobilization of 
HRC trimer peptide. As shown in Fig. 12 and Table 4, 
except for HRC7, HRC, and HRC analogs could bind the 
HRC trimer chip, which is consistent with native-PAGE 


analysis. 
Based on the results in Table 4 obtained from CD, the 
complex stability order of HRC/HRC analogs is 


HRC4> HRC3> HRC7 ~ HRC1> HRC5>HRC2>HRC= 
HRC6. All analogs showed interaction with HRC by native gel 
electrophoresis and SPR analysis, except for HRC7. HRC6 
showed very weak interaction with HRC since the complex 
has similar helical structure and stability (T,,=5°C) to 
native HRC. Therefore, the interaction of the analogs with 
native HRC correlates with stability, except for HRC7. 


4. Discussion 


Formation of the six-helix bundle is a critical conforma- 
tional change of the S protein during fusion of the SARS- 
CoV and target cell membranes. Clearly therefore, the 
functional domains of the S protein involved in membrane 
fusion represent attractive targets for the discovery of viral 
entry inhibitors of SARS-CoV. A compound that blocks 
the S protein six-helix bundle formation may be a potent 
SARS-CoV fusion inhibitor. Therefore, determination of 
S protein six-helix bundle formation is important for 
studying the mechanism by which SARS-CoV induces 
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Fig. 12. Surface plasmon resonance analysis of the interaction of HRC3 with HRC trimer chip. HRC trimer was covalently coupled to the Biacore CM5 
chip surface. Binding experiments were performed at 25 °C in PBS buffer, pH 7.4, with 0.005% Tween 20. The peptides to be tested w ere introduced at four 
serial concentrations (5, 10, 20, and 50 1M) in running buffer in 80 pl injection volumes and at a flow rate of 20 l/min, followed by 6 min of dissociation. 
HRC3 is shown in the schematic as a monomer-trimer equilibrium in solution when flowing over the immobilized HRC trimer on the chip surface. 
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membrane fusion and for identification of SARS-CoV 
fusion effective inhibitors targeting S protein. In this bio- 
physical study, we made substitutions to HRC and studied 
their effect on o-helical structure and stability in the 
absence and presence of the native sequences of HRN and 
HRC. 

In the crystallographic structures of the six-helix bundle 
conformation (Duquerroy et al., 2005; Supekar et al., 2004; 
Xu etal., 2004a), the 49 residues (902-950) in the HRN 
region and the 36 residues (1150-1185) in the HRC region 
were identified to be the major binding regions for their 
interactions. Interestingly, the 36-residue length of HRC 
matches the 49 residues of HRN in a trimeric coiled-coil 
conformation. This is accomplished by residues 1161-1175 
in HRC forming a 5-turn o-helix with residues 1150-1160 
at the N-terminus and residues 1176-1185 at the C-termi- 
nus of HRC forming two extended conformations. This 
result suggests the N- and C- terminal residues of HRC 
must be able to readily change conformation from an o- 
helical conformation in a trimeric coiled-coil to the 
extended conformation when bound to HRN in the six- 
helix bundle conformation. 

Substitutions in the region of HRC could have one of 
three effects on its interaction with HRN: first, substitution 
of a particular native residue could affect the interaction of 
that side-chain with HRN, thus decreasing the stability of 
the HRN/HRC complex; second, substitution could stabi- 
lize the o-helical conformation in HRC which in turn 
would reduce the ability of the N- and/or C-terminal 
regions of HRC to readily change conformation from o- 
helical to extended on interacting with HRN; third, if the 
monomeric form of HRC is the interacting species with 
HRN, then substitutions that stabilize the oligomeric struc- 
ture would shift the monomer—oligomer equilibrium in the 
favor of the oligomer and at some point no monomer 
would be available to interact with HRN. 

In our analysis of the interaction between HRC and 
HRN, we used the X-ray crystallographic structure of the 
six-helix bundle conformation of Supekar et al. (2004) (pdb 
accession number is | BEZ) and for the interactions in HRC 
we used the solution NMR structure of residues 1141-1193 
by Hakansson-McReynolds et al. (2006) (pdb accession 
number 2FXP). Our choice of the HRC region for inhibitor 
design (HRC 1150-1185) was initially based on our bio- 
physical studies (Tripet etal. 2004) before any three- 
dimensional structural information was available. This 
choice was confirmed by the NMR solution structure which 
showed that the trimeric coiled-coil was stabilized by a 
hydrophobic core involving hydrophobic residues in the a 
and d positions from Ile 115la to Leu 1182d (Fig. 1). The 
NMR structure also showed the importance of Leu 1185g 
in hydrophobic packing between the three polypeptide 
chains. The importance of hydrophobes in positions e and g 
of coiled-coils for increasing coiled-coil stability has been 
previously documented by biophysical studies and the X- 
ray crystallographic structure of a cortexillin I/GCN4 
hybrid coiled-coil (Lee et al., 2003). 


In the six-helix bundle conformation the interaction 
between HRC and HRN can be divided into three regions 
based on HRC (Fig. 1). Region 1 involves residues 1150- 
1160 of HRC which are helical in the prefusion state (NMR 
solution structure of the HRC trimer) and extended in the 
6-helix bundle conformation. That is, in this region hydro- 
phobes Ile 1151a, Ile 1154d, and Val 1158a are in the hydro- 
phobic core in the pre-fusion state and hydrophobe Val 
11595 is exposed on the outside of the trimeric coiled-coil 
of HRC. In contrast, in the six-helix bundle conformation 
Val 11596 is now buried along with Ie 1151a and Ile 1154d 
on interaction with HRN and Val 1158qa is now exposed. 
That is, there has been a switch of hydrophobes between 
the two states involving Val 1158a and Val 11595. Region 2 
involves residues 1161-1175 which are helical in both the 
prefusion state (HRC trimeric state) and the six-helix bun- 
dle conformation and the same four hydrophobes (Ile 
1161d, Ile 1165a, Leu 1168d, and Leu 1175d) are stabilizing 
this helical structure in both states (interacting with the 
other chains of the HRC trimer or interacting with the 
HRN trimer). In this region there is also one hydrophobe 
not in the a and d positions, Val 1171g, which is interacting 
to stabilize the hydrophobic core in the HRC trimer and is 
also involved in the interaction with HRN. This residue is 
conserved as a hydrophobe in all coronaviruses as Val, Leu 
or Ala (Duquerroy et al., 2005). Region 3 involves residues 
1176-1185 of HRC which are helical in the HRC trimer 
and extended in the six-helix bundle conformation. That is, 
hydrophobes Leu 1179a, Leu 1182d, and Leu 1185g are in 
the hydrophobic core of the HRC trimer and Ile1180b is 
exposed on the outside of the trimeric coiled-coil of HRC. 
In contrast, in the six-helix bundle conformation Ile 1180b 
is now buried along with Leu 1182d on interaction with 
HRN, and Leu 1179a is now exposed (Leu 1185g is not 
present in the crystal structure of the six-helix bundle). That 
is, there has been a switch of hydrophobes between the two 
states involving residues Leu 1179a and Ile 1180b. It is 
interesting that the two hydrophobes, Val 11594, and Ile 
11804 that are surface exposed in the HRC trimer and bur- 
ied on interaction with HRN in the six-helix bundle confor- 
mation are conserved among all coronaviruses as a 
hydrophobe (Val or Leu at position 1159 and a B-branched 
hydrophobe Val or Ile at position 1180), indicating their 
importance to presumably stabilize the six-helix bundle 
conformation (Duquerroy et al., 2005). 

The question to be addressed is whether we can explain 
the behaviors of our HRC analogs in light of the structural 
information. A HRC inhibitor could inhibit six-helix bun- 
dle formation in several ways: by binding to native HRC, 
by binding to native HRN or by interacting with both 
regions. The requirement for the complete HRC sequence 
1150-1185 for binding to HRN was demonstrated by our 
deletion analogs (Table 3). The importance of the N-termi- 
nal region of HRC to interact with HRN was demonstrated 
by the 25-residue peptide (HRC 25) which had 11 residues 
(1150-1160) deleted from HRC was unable to interact with 
HRN. Adding back 3 additional residues (VVN, residues 
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1158-1160) to HRC 25 produced HRC 28 which bound to 
HRN witha T,, value for the HRC/HRN complex of 58 °C, 
compared to the stability of HRN-HRC complex of 85°C 
(Table 2). This result illustrates that the N-terminal seven 
residues (1150-1157) are extremely important for HRC 
binding to HRN. In particular, based on the X-ray struc- 
ture of the six-helix bundle Ile 1151 or Ile 1154 interact with 
HRN in the hydrophobic core as discussed above. Simi- 
larly, truncation of the N-terminal residues 1150-1160 from 
the HRCS sequence (Fig. 1), the best binding analog to 
HRN (Table 2) also showed no interaction with HRN. This 
result is also supported by Duquerroy et al. (2005), who 
reported that the extended regions of HRC consolidate 
their interaction with HRN by a string of asparagine and 
glutamine side chains from the HRN coiled-coil, which 
hydrogen bond to the HRC main chain. This extensive 
hydrogen-bonding network zips the HRC chain along the 
HRN interhelical grooves. 

Four HRC analogs (HRCI, HRC2, HRC3, and HRC4, 
Fig. 1) were designed to increase o-helical structure and sta- 
bility of the o-helix and resulting coiled-coil by either sub- 
stitution of residues in positions b, c, and f of the heptad 
repeat by Ala (outer positions on the coiled-coil) to increase 
a-helical propensity or substitution of Ala by Ile in position 
a to increase hydrophobicity in the hydrophobic core. The 
order of _ stability of these analogs was 
HRC2< HRCI <HRC3<HRC4 (Table 1). Interestingly, 
this identical order was obtained for the stability of 1:1 
mixtures of these analogs with HRC (Table 4). These 
results suggest that increasing the stability of the HRC ana- 
logs increases their ability to interact with the native HRC 
sequence. The most stable analog of this series of analogs 
was HRC4 with a T., of 74°C compared to HRC with a T,, 
of 37°C at a peptide concentration of 70uM (Table 1). 
Thus, this dramatic increase in stability does not prevent 
interaction with HRC and suggests that the stability of this 
HRC analog is not so high as to prevent interaction with 
the native HRC sequence whether the interaction is based 
on the oligomeric coiled-coil or monomer of the analog 
interacting with the native HRC sequence. 

Interestingly, HRC and three HRC analogs that are the 
least stable oligomeric coiled-coils (T,,, values of 37, 48, 41, 
and 32°C for HRC, HRC2, HRCS5, and HRC6, respec- 
tively, Table 1) formed the most stable complexes with 
HRN (7,, values of 85, 86, 88, and 83°C for HRN com- 
plexes involving HRC, HRC2, HRCS, and HRC6, respec- 
tively, Table 2). These results suggest that you can increase 
the stability of the oligomeric form of HRC with little effect 
on the formation and stability of the HRN/HRC complex 
(11°C increase in T,, from 37°C for HRC to 48°C for 
HRC2). On the other hand, the substitutions in HRC which 
stabilized the oligomeric coiled-coil structure the most 
(T,, values of 57, 66, and 74°C for HRC1, HRC3, and 
HRC4, respectively compared to HRC with a T,,, of 37°C, 
Table 1) formed the least stable complexes with HRN 
(T,, values of 74, 75, and 59°C, respectively, compared to 
HRC/HRN complex with a T,, of 85°C, Table 2). 


The substitutions that increase hydrophobicity in the 
hydrophobic core (Ala 1172 Ile) increased stability of the 
trimer (HRCI and HRC3, Table 1) as expected. These two 
analogs also bind to native HRC (Table 4). However, this 
Ala to Ile substitution was disruptive to the stability of the 
HRN/HRCI1 or HRN/HRC3 complex, compared to the 
native HRN/HRC complex (Table 2). Obviously, Ile is 
more ideal for stabilizing the trimer than Ala but Ala at this 
position is more ideal than Ile for stabilizing the HRN 
interaction. This suggests that though this region 2 (Fig. 1) 
is helical in the trimer and helical on interacting with HRN, 
the interactions between the helix and HRN are different 
from the helix-helix interactions in the trimer. 

HRCS has a i to i+4 covalent lactam bridge between 
Glu and Lys (Fig. 1) which did provide an increase in sta- 
bility of the trimeric coiled-coil (compare the control pep- 
tide HRC6 with a salt bridge at the same positions with a 
T,, of 32°C to HRCS5 with a T,, of 41°C). However, the 
substitutions of Asp and Glu in HRC to Glu and Lys in 
HRC6 suggest that Asp and Glu must be important for the 
stability of the trimeric coiled-coil of HRC since HRC6 was 
less stable than the native HRC sequence by 5°C (Table 1). 
This clearly demonstrates that the covalent constraint does 
stabilize the o-helix and the interaction of HRCS5S with 
HRC but that the location of the constraint should be such 
as to not disrupt the coiled-coil in the control peptide with 
the non-covalent ion-pair (HRC6). Thus, we probably did 
not achieve maximum benefit of this covalent constraint. 
Since HRCS stabilized the complex with HRN more than 
any other analog (Table 4), HRCS5 should be classified as a 
potential HRN-specific or combined HRN and HRC 
inhibitor. Thus, clearly stabilizing the a-helical region in the 
six-helix bundle state is advantageous. 

In fact, HRC4 (8 Ala substitutions to increase o-helical 
propensity, Fig.1) formed the most stable oligomeric 
coiled-coil of all the HRC analogs (T,, of 74°C compared 
to 37°C for native HRC, Table 1) and the least stable com- 
plex with HRN (T,, of 59 or 26°C below the stability of the 
native HRN/HRC complex, Table 2). On the other hand, 
HRC4 formed the most stable complex with native HRC 
increasing the T., by 50°C (Table 4). These results can be 
explained based on structural information. Stabilizing the 
helix analog along the length of the trimeric coiled-coil pro- 
vides an excellent potential inhibitor of the six-helix bundle 
state by complexing with native HRC (a potential HRC- 
specific inhibitor), but stabilizing region 1 and 3 (Fig. 1) in 
the helical conformation prevents the conformational 
change of these regions to go to the extended conformation 
in the six-helix bundle conformation. Also, we have substi- 
tuted the two conserved residues Val 1159 and Ile 1180 with 
Ala and these residues are critical for interaction with 
HRN. Thus, HRC4 would not exert its inhibitory effect by 
binding to HRN. 

HRC7, with 10 substitutions of Glu and Lys into the 
HRC sequence at positions 5, c, f, and g to create a series of 
ito i+3 and i to i+4 intrachain electrostatic interactions 
(Fig. 1), did stabilize the o-helical structure of the trimer by 
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24°C compared to native HRC (Table 1). CD results sug- 
gested an interaction of HRC7 and HRC (Fig. 10) but no 
interaction was observed by Native-PAGE (Fig. 11) or SPR 
(Table 4). Though these salt-bridges stabilize the coiled-coil 
structure, at least some of these side-chains in this analog 
must be critical for disrupting interactions with native 
HRC. Thus, the interaction with HRC requires more than 
just the packing of the a and d positions of the hydrophobic 
core which are identical in HRC7 and native HRC. In addi- 
tion, HRC7 showed absolutely no binding to HRN. This 
can be explained by the fact that the 3 conserved hydro- 
phobes in all coronaviruses (positions 1159, 1171, and 
1178) that are important for binding to HRN have been 
substituted by charged residues which would disrupt and 
destabilize the interaction with HRN. Recently it was also 
reported that position 1160 is a conserved Asn/Asp residue 
in coronaviruses and from the structure of the six-helix 
bundle the side-chain hydrogen bonds to the polypeptide 
backbone of HRC and Glu 1164 is conserved for interac- 
tion with HRN via a salt bridge (Duquerroy et al., 2005). 
Therefore, the substitution of these five positions would be 
expected to abolish binding to HRN. 

In summary, the above results show that increasing a- 
helical propensity of side-chains, increasing hydrophobicity 
in the hydrophobic core, increasing intrachain ionic inter- 
actions or addition of a covalent constraint that increase 
coiled-coil stability over the native sequence can be benefi- 
cial to the design of a potential peptide inhibitors that are 
HRC-specific (HRCI, HRC3, and HRC4) or inhibitors 
that can bind both HRC and HRN simultaneously (HRC2 
and HRC5). Biological experiments will be carried out to 
verify the inhibitory activity of these potential peptide 
inhibitors and whether HRN or HRC is the better target. 
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